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ABSTRACT 


A  study  of  two  oceanic  extra  tropical  cyclones,  generated 
by  the  Navy  Operational  Atmospheric  Prediction  System 
(NOGAPS)  model  (see  Sandgathe,  1981),  ravaals  a  simultaneous 
development  of  the  upper-  and  lower-tropospheric  features 
while  these  features  maintain  their  positions  relative  to 
each  other.  These  cyclones  also  develop  extreme  shear  in 
most  synoptic  fields  after  120  hoars  of  numerical 
simulation. 

The  mass  budgets  of  these  cyclones  show  a  concentration 
of  inward  lateral  transport  (oonvergencs)  in  the  lowest 
model  layer  and  outward  lateral  transport  (divergenca)  in  a 
layer  centered  near  300  mb .  Time  sections  of  lateral  trans- 
port, vertical  velocity  and  mass  tendency  reveal  that  these 
cyclones  develop  in  two  phases — an  explosive  cyclogenesis 
phase,  and  a  quasi-steady  state  phase.  These  time  sections 
also  reveal  a  12-hour  cyclic  pattern  embedded  within  the  two 
phases  of  cyclone  evolution.  This  12-hour  cyclic  pattern  is 
not    evident  in   the   map    sequences    of  the   surface   pressure. 
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I-       121X13 SUCTION 

For  many  years  maritime  cyclones  have  been  of  interest 
to  military  and  civilian  communities  from  both  a  scientific 
and    a     practical   point      of    view.  Oceanic   lows      frequently 

develop  rapidly  and  produce  meteorological  phenomena  which 
can  adversely  impact  civilian  and  military  operations  in 
oceanic  and  coastal   regions.  Accurate   and    dependable  pre- 

diction of  the  development  and  movement  of  oceanic  cyclones 
is    required   for   reliable   maritime   forecasts. 

Numerous  studies  of  middle-latitude  cyclones  have 
focused     on        continental     cyclones      (e.g.  Smith      (1973)  , 

Petterssen  and  Smebye  (1971),  Johnson  and  Downey  (1976), 
Chen  and  Bosart  (1977),  and  others).  This  is  primarily  due 
to  the  abundant  observational  data  which  is  readily  avail- 
able from  continental  regions.  In  contrast,  relatively  few 
studies  have  investigated  maritime  cyclones  because  of  the 
limited  observational  data  from  oceanic  regions.  As  a 
result,  the  factors  influencing  maritime  cyclogenesis  are 
not    fully   understood. 

Several  studies  of  oceanic  cyclones  indicate  that  sig- 
nificant differences  exist  in  the  physics  involved  in  conti- 
nental and    oceanic   cyclogenesis.         Petterssen,      Bradbury   and 


12 


Pedersen  (1962)  compared  the  development  of  cyclones  over 
the  western  North  Atlantic  to  that  of  cyclones  over  the 
eastern  United  States.  This  observational  study  indicates 
that  the  initial  development  of  oceanic  cyclones  appeared  to 
result  primarily  from  low- lev  el  thermal  advection  while  con- 
tinental cyclones  usually  developed  as  a  result  of  an 
upper- level  vorticity  maximum  approaching  a  low-level  baro- 
clinic  zone.  A  study  by  Nitta  and  Yamamoto  (1974)  also  sug- 
gests that  maritime  cyclones  frequently  develop  in  the 
vicinity  of  an  ocean  baroclinic  zone  ia  conjunction  with 
nearly  zonal  upper- level  flow.  Sanders  and  Gyakum  (1980) 
found  that  explosive  cyclogenesis  is  predominantly  a  mari- 
time, cold-season  event  and  suggested  that  the  roles  played 
by  cumulus  convection  and  air-sea  fluxes  in  cyclogenesis  are 
not  adequately  understood.  A  numerical  study  by  Sandgathe 
(198  1)  indicates  that  air-sea  fluxes  play  a  significant  role 
in  the  evolution  of  oceanic  cyclones.  These  studies  and 
others  suggest  that  the  air-sea  fluxes  of  latent  and  sensi- 
ble heat  play  a  significant  role  in  oceanic  storm 
development. 
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Current  operational  numerical  models  have  been  shown  to 
be  seriously  deficient  in  predicting  maritime  cyclone 
deepening  rates.  The  Sanders  and  Gyakum  (1980)  study  indi- 
cated that  the  National  Meteorological  Canter  (NMC)  primi- 
tive-equation models  predicted  less  than  one  half  of  the 
deepening  rates  observed  in  explosive  oceanic  cyclogenesis. 
The  failure  of  numerical  models  to  predict  explosive  cyclo- 
genesis is  not  entirely  the  result  of  the  limited  observa- 
tional data  from  oceanic  regions.  Part  of  the  problem  lies 
in  modeling  vertical  exchanges  in  the  atmosphere.  The  upper 
troposphere  is  dominated  by  wave  structures,  while  vortex 
structures  dominate  in  the  lower  troposphere.  Most  cases  of 
upper  tropospheric  wave  amplification  are  accompanied  by  the 
development  of  anticyclones  and  cyclones  in  the  lower  tropo- 
sphere. This  suggests  that  there  exists  a  strong  interacw 
tion  between  these  two  regimes.  The  smaller  size  of  oceanic 
cyclones  presents  another  problem.  The  space  scale  of  upper 
tropospheric  waves  is  typically  two  to  four  times  the  space 
scale  of  oceanic  cyclones.  Us  a  result,  it  appears  that 
numerical  models  adequately  simulate  the  large-scale 
processes  involved  in  synoptic-scale  systems,  but  do  not 
capture        the        smaller-scale      features        embedded  within 

synoptic-scale   weather   systems. 
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Due  to  the  sparse  observational  network  in  oceanic 
regions,  few  quantitative  studies  of  maritime  cyclogenesis 
have  been  completed.  To  fill  this  void,  studies  of  numeri- 
cally simulated  cyclone  development  are  needed-  Numerical 
models,  having  an  internally  consistent  "data"  base,  elimi- 
nate errors  resulting  from  incomplete  or  inaccurate  observa- 
tions from  data-sparse  regions  such  as  oceans.  In  cases 
where  there  is  a  favorable  comparison  between  numerically 
simulated  cyclogenesis  and  real  atmosphere  cyclogenesis,  the 
model  may  provide  a  greater  un derstanding  of  the  physics 
involved  in  cyclone  development  than  would  the  incomplete 
observations  of  the  real  atmosphere  storm.  A  disagreement 
between  the  model  and  the  real  atmosphere  would  not  neces- 
sarily prove  to  be  of  little  or  no  value.  In  cases  where 
the  model  and  the  real  atmosphere  do  not  agree,  the  diagnos- 
tics of  the  real  atmosphere  may  well  isolate  the  reasons  for 
the    disagreement   and   thus   lead  to   improvement   in   the   model. 

A  basic  feature  of  cyclogenesis  is  an  inward  mass  trans- 
port (convergence)  in  the  lower  troposphere  and  and  outward 
mass  transport  (divergence)  in  the  upper  troposphere.  Asso- 
ciated with  these  horizontal  mass  transports  is  an  upward 
vertical    mass      transport   within      the   cyclone.  Quantifying 
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this  mass  circulation  can  provide  important  information  about 
the    large-scale    forcing   within   a    cyclone. 

Cyclones  are  vortices  embedded  in  a  planetary-scale 
pattern.  Consequently,  a  study  of  the  basic  properties  of 
cyclones  should  take  advantage  of  their  vortex  nature.  One 
such  approach  is  the  investigation  of  the  budgets  of  the 
fundamental  quantities  such  as  mass  and  circulation  within  a 
volume  moving      with  the    storm.  Quasi-Lagrangian   transport 

theory,  outlined  by  Johnson  and  Downey  (1975  a  and  b) ,  is 
particularly   suited   for   the  study  of   cyclones. 

In  this  thesis  the  mass  budgets  of  nunerically  simulated 
oceanic  cyclones  are  investigated  to  ascertain  the  basic 
aspects   of   cyclone  development   and    maintenance.  A   quanti- 

tative study  of  the  mass  circulation  helps  describe  the 
intensity  and  physical  structure  of  the  cyclone  and  isolate 
vertical    variations   in   the   cyclone  circulation. 

This  study  follows  the  Ph.  D.  thesis  investigation  by 
Sandgathe  (1981).  The  budget  studies  were  applied  to  two 
storms  generated  by  a  version  of  the  CJCLA  general  circula- 
tion model,  which  includes  a  sophisticated  parameterization 
scheme  for  planetary  boundary  layer  (PBL)  and  free  atmo- 
sphere interactions-        The   Sandgathe   study      demonstrated  the 
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role  of  air-sea  fluxes  of  heat  and  moisture  in  initial 
cyclogenesis.  He     showed       that     cyclone      intensification 

resulted  when  sensible  heat  fluxes  were  removed.  In  con- 
trast, the  lack  of  moisture  fluxes  resulted  in  cyclone 
decay. 

In  this  thesis,        two    cyclones   were  studied — one      in  the 
southern     (spring)      hemisphere   and  one   in   the   northern     (fall) 
hemisphere.  Both        storms      are     characterized        by      rapid 

development  during  the  period  from  96  to  126  hours  of  simu- 
lated cyclogenesis.  In  each  case,  the  central  surface  pres- 
sure of  these  storms  decreases  at  a  rate  of  nearly  one  mb 
per  hour  during  this  30  hour  period.  Thus,  these  storms  are 
good  examples  of  the  explosive  deepening  which  is  not 
adequately    handled   by   current   operational   dynamic  models. 

This  study  is  one  in  a  series  of  experiments  designed  to 
investigate  the  physical  processes  involved  in  extratropical 
cyclone  development.  The  primary  objectives  of  this  study 
are: 

•  to  study  the  synoptic  features  of  two  model-generated 
cyclones  and   compare   them   to  observed   cyclones; 

•  to  isolate  the  lateral  and  vertical  branches  of  the 
mass  circulation  during  different  stages  in  the 
evolution   of   a    model- genera  ted   maritime    cyclone;    and, 

•  to  compare  vertical  velocity  values  derived  through 
quasi-Lagrangian  diagnostics  with  the  vertical  velocity 
values  from  the  models  internally  consistent  "data" 
set . 
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The  atmospheric  model  used  for  this  study  is  briefly 
described  in  Chapter  II.  The  rssults  of  this  study  begin  in 
Chapter  III  with  an  objective  review  of  the  synoptic  fea- 
tures of  the  two  model  storms.  A  description  of  the  guasi- 
Lagrangian  diagnostics  technigue  and  the  eguations  used  in 
the  mass  and  circulation  budgets  are  prssented  in  Chapter 
IV.  The  generalized  budget  eguation  and  the  budget  volume 
coordinate  system  are  described  in  Appendix  B.  Chap+er  7 
contains  the  analysis  of  the  mass  budget  results.  Conclu- 
sions and  recommendations  for  further  studies  are  presented 
in  Chapters  VI  and  VII  respectively. 
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II.       ATMOSPHERIC    M3DEL    DESCRIPTION 

A.        MODEL    STRUCTURE 

The  Sandgathe  experiments  used  the  Navy  Operational 
Atmospharic  Prediction  System  (NOGAPS)  madel  which  is  being 
developed  by  the  Naval  Environmental  Prediction  Research 
Facility    (NEPRF) .  The   NOGAPS    model   is      a    state-of-the-art 

version  of  the  UCLA  general  circulation  model.  A  brief 
overview  of  the  NOGAPS  model  is  presented  in  this  chapter. 
For  a  more  detailed  description  of  the  NOGAPS  model  the 
reader  is  directed  to  Sandgatha  (1981).  The  UCLA  model  is 
described    in   detail   by   Arakawa   and  Lamb    (1977). 

A  fine  resolution  grid  (2  degress  latitude  by  2.5 
degrees  longitude)  within  a  5  0  degree  longitude  global 
sector  is  used  in  this  model.  This  grid  resolution,  which 
eguates  to  a  222  km  by  196  km  grid  distance  at  45  degrees 
latitude,    is   sufficient    to   resolve   shorter   cyclone  waves. 

The  prognostic  variables  of  the  model  include  the  sur- 
face pressure  minus  50  mb,  tt  ,  the  horizontal  wind  velocity, 
V,  temperature,  T,  and  specific  humidity,  q.  Several  other 
prognostic  variables  (to  be  discussed  latar)  are  associated 
with   the    PBL. 
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The      model      uses      a      sigma      vertical      coordinate      system 
defined   as: 


IT 

where: 

P   =    pressure 
IT  =    Ps    -    PT 

Ps    =   surface   pressure 

PT  =  pressure  at  the  top  of  the  model  atmosphere  (50  mb)  . 
Thus  the  top  and  bottomy  moundaries  of  the  model  ars  a  -  0 
and  (7=1.  The  NOGAPS  model  consists  of  six  vertical  layers 
(Fig.  1)  .  In  this  depiction,  an  arbitrary  surface  pressure 
of  1000  mb  has  been  assumed.  Variables  carried  at  the  layer 
mid-levels  are  T,  V,  and  g.  The  vertical  velocity,  <7,  is 
defined   at   the  layer   interfaces. 

The  model  employes  a  staggsred  grid  (Arakawa  scheme  C) 
in  the  horizontal  (Fig.  2).  The  variables  defined  at  the 
canter  points  are  T,  g,  and  tto .  The  msridional  wind  compo- 
nent, v,  is  carried  one  degree  north  and  south  of  the  center 
point  with  the  zonal  wind  component,  u,  carried  1.25  degrees 
east  and  west  of  the  center  point.  For  diabatic  and  fric- 
tional  computations,  the  u  and  v  components  are  averaged  to 
the  center  point.  A  smoothing  of  the  msridional  component 
adjacent    to   the    poles   is   implemented    in    theis    model. 
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Figure   1.      Vertical   distribution    of   model   large-seals 

prognostic   variables.       Pressurs    values   of   sigma 
levels    vary    with  tha    surface   prassure.      A   surface 
value   of    1000    ab  is   assumed   in    this    fiaure. 
(After   Sandgathe,    198  1) 
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The  area  used  for  this  study,  50  degrees  wide  extending 
from  pole  to  pole,  is  an  open-ocean  (i.e.  entirely  over 
water  with  no  sea  ice)  domain  which  is  fixed  in  time  with 
two  zones  (centered  at  U5N  and  35S)  of  maximum  north-south 
gradients  in  sea-surface  temperature.  Continuous  over-ocean 
propagation  within  the  model  is  guaranteed  by  imposing 
cyclic  boundary  conditions  in  the  east-west  direction. 
There  are  91  and  24  grid  points  in  the  north-south  and 
east-west  directions  respectively.  For  this  thesis,  zero 
degrees  longitude  has  arbitrarily  been  assigned  to  the 
western    boundary. 

The  PBL  paramererization  scheme  of  the  NOGAPS  model, 
which  follows  Deardorff  (1972)  and  Randall  (1976),  provides 
for  interactions  between  the  PBL  and  cumulus  and/or  stratus 
cloud  ensembles  at  each  grid  point.  Momentum,  heat  and 
moisture  fluxes  at  the  surface  are  diagnos tically  determined 
using  the  values  of  T,  7,  and  q  from  the  dynamic  portion  of 
the  model  and  a  bulk  Richardson  number  which  is  based  no 
sea-surface     temperature.  The        cumulus      parameterization 

scheme  follows  Arakawa  and  Schubert  (1974}  and  Lord  (1978). 
Effects  of  water  vapor,  predicted  cloud  distributions,  and 
large-scale  precipitation  processes  are  included  in  the 
radiational   heating   package   which   follows   Katayama    (1972). 
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Figure  2.      Horizontal  distribution   of   modal    large-scale 

Drognostic   variables.       This   distribution   follows 
Arakawa   scheme   C.       (After   Sandgathe,    1981) 


B.        INITIAL    CONDITIONS 

The  Sandgathe  study  utilized  initial  conditions  which 
are  representative  of  open-ocean  fall  and  spring  regimes. 
The  initial  conditions  for  the  model  atmosphere  follow 
Simmons      and      Hoskins       (1977,    1978).  Wind        profiles      for 

the      fall      hemisphere      were     determined   using      a      polynomial 
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approximation  of  the  January  mean  zonal  wind  profile  at  30N 
from  Oort  and  Rasmussen  (1971).  The  fall  hemisphere  winds 
are  multiplied  by  a  factor  of  1.5  at  each  level  to  provide 
the    spring   hemisphere   wind   profiles. 

The  model  start  time  was  1  Dctober.  Thus,  with  respect 
to  solar  fluxes,  the  northern  and  southern  hemispheres  are 
the    fall      and      spring     regimes      respectivsly.  Sea-surface 

temperatures  and  atmospheric  jets  are  also  representative  of 
these  seasons.  Zonal    jets   of    35      m   sec-1    and   52      m    sec-1, 

centered  at  U5N  and  35S  respectively,  represent  possible 
conditions    for  these     regimes.  As      shown      by      Simmons   and 

Hoskins  (1977,  1978),  the  existance  of  thase  jets  guarantees 
the    growth   of  baroclinically    unstable    waves. 

Sandgathe  used  the  conditions  listed  above  as  initial 
conditions  in  his  experiments  on  the  evolution  of  extratrop- 
ical  cyclones.  A  small  amplitude  sinusoidal  perturbation  in 
wave  number  six  was  added  to  the  meridional  wind  component 
at  each  level.  The  maximum  value  of  the  perturbation  was  .5 
m  sec-1  at  U5N/S.  The  perturbation  at  the  equator  and  the 
poles  was  zero.  The  associated  perturbation  in  the  surface 
pressure    field  was   approximately    one    millibar. 
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III.       SYN0PTI3    DISCUSSION 

A.        OVERVIEW 

An  initial  objective  of  this  thesis  is  a  comparison  of 
the  synoptic  features  of  numerically  simulated  and  observed 
oceanic  cyclones.  This  discussion  will  focus  on  the  signif- 
icant features  in  the  model  "data"  fields.  Available  "data" 
consists  of  model  output  at  three-hour  intervals  from  hour 
96  to  hour  168  inclusive.  Except  for  the  1000  mb  heights 
and  surface  pressures,  all  fields  in  this  discussion  are 
generated  from  mid-layer  model  data.  With  the  exception  of 
the  wind  fields,  a  three  point  zonal  smoothing  technique  was 
applied  to  all  data  fields  raprasented  in  this  chaptar.  The 
u  and  v  wind  components  were  averaged  to  the  height  grid 
points  for   the   wind   field   representations. 

The  synoptic  summary  of  tha  two  model  cyclones  begins 
with  a  presentation  of  the  significant  features  in  the  lower 
troposphere  using  surface  pressure  and  temperature  fields, 
the  1000  mb  height  fields  and  the  1000  mb  to  mid-layer  U 
(approximately  1000-500  mb)  thickness  fields  (Figs.  4  to 
11) .      This    is      followed      by      a      summary      of    mid-tropospheric 


25 


features  (aid-layer  4,  approximately  530  mb)  of  the  height, 
temperature,  wind,  and  vorticity  fields  (Figs.  12  to  19). 
Significant  features  in  the  mid-layer  3  (approximately  300 
mb)  height,  wind,  and  temperature  fields  are  described  in 
the  fourth  section  (Figs.  20  to  27).  This  chapter  concludes 
with  a  brief  comparison  of  the  development  of  the  model 
cyclcnes    and   observed  cyclones. 

Table  I  contains  a  list  of  synoptic  fields  relevent  to 
this  summary.  Fields  valid  at  96,  120,  144,  and  168  hours 
after  the  initial  time  are  presented.  With  respect  to  the 
synoptic  summary  and  the  figuras  presented  in  this  thesis, 
the  indicated  times  represent  the  number  of  hours  of  numer- 
ical simulation.  The  storm  tracks  of  the  two  numerically 
•simulated  cyclones  used  for  the  budgst  studies  of  this 
thesis  are   shown   in   Fig.    3. 

The  figures  presented  in  this  chapter  depict  two 
adjacent  60  degree  sectors  based  on  the  cyclic  continuity 
assumption  of  the  model.  Six  and  nine  dsgree  budget  volume 
circles,  centered  on  the  surface  prassure  minimum,  are 
depicted    on     each   chart.  Radii   six      and    nine      were   chosen 

because  they  include  the  more  intenss  portion  of  the 
cyclcnes    during      the   early      and    latter      stages   of     the    study 
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respectively.  In  following  discussions,  inner  radius  and 
outer  radius  refer  to  the  six  and  nine  degree  budget  radii 
respectively.  Both  the  northern  hemisphere  (NH)  and  south- 
ern hemisphere  (SH)  cyclones  are  discussed,  relative  to  each 
figure,  to  point  out  differences  and  similarities  in  their 
development. 


TABLE    I 
DATA    FIELDS    DEPICTED    IN    THE    SYNOPTIC    SUMMARY 


SIGMA  APPROXIMATE 

LEVEL  PRESSURE 

(OR    LAYER)  LEVEL(S)  FIELDS 


level    6  surface  pressure    (mb)  / 

temperature 

(1000    mb)/  1000/  1000    -   500    mb 

mid-layer   4  500   mb  thickness/ 

1000    mb   height 

level    6  surface  pressure/ 

temperature 

mid-layer   4  500    mb  vorticity/ 

height 

mid-layer   4  500   mb  temperature/ 

height 

mid-layer   3  300    mb  temperature/ 

height 

aid-layer   3  300    mb  w^nd   vectors/ 

'  isotachs 
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Figure  3. 


Storm   tracks   of   two  model    cyclones.      Indicated 
locations   are   for   six    hoar   intervals    from   hour    96 
through   hour    168.      Circle   radius    is   6°    latitude, 
a.    Northern   hemisphere.      b.    Southern   hemisphere. 


B.        LOWER-TROPOSPHERE    FEATURES 

Surface  pressure  and  temperature  fields  valid  at  hour  96 
are  shown  in  Fig.  4.  The  100  0  mb  height  and  1000  mb  to 
mid-layer    4    (approximately    1000    -     500   mb)       thickness    fields 
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for    the      same  time  are      presented    in   Fig.    5.  Wealc   thermal 

ridges,      at  32    (92)      and   57    (117)      degrees    longitude   respec- 
tively,   are  apparent    in   the  NH.         A   1011    ib   low    (cyclone   N1) 
is    located   in   the  thermal   ridge    near    32    (92)       degrees    longi- 
tude.     A    somewhat   more    intense  low    (cyclone   N2)  ,    with    a  cen- 
tral  surface   pressure      of    1007   mb,      is   in      the   thermal    ridge 
naar   57       (117)       degrees   longitude.        Initially,        cyclone    N2 
appears  to   be  the   main      feature;       however,       within   2<*    hours, 
cyclone   N1      becomes  the   dominant    low      in   the    NH      sector.         A 
third  area   of      low   pressure  approximately    1000      nm   northeast 
of    cyclone  N1    develops   very  little,    but   does    persist   through 
hour    156.        This   low    pressure   area   will   be    discussed   briefly 
later  in    this   section.      In   the   SH,   a    1003   ib   low    (cyclone    S) 
is    located    in   a    thermal    ridge   oriented   south-southwest     (SSW) 
to    north-northeast    (NNE)    and   centered  near    40    degrees    longi- 
tude.     Table   II   contains   the   central    surface    pressure    values 
of    the      three   cyclones   discussed      in   this      synoptic   summary. 
Cyclones        N1,        N2,      and      S      all     display        typical      lower- 
tropospheric   thermal   advection   patterns   with      warm   advection 
ahead   of,    and  cold   advection   behind      and      eguatorwardof   the 
developing   cyclones. 
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TABLE    II 
CENTRAL    SURFACE    PRESSURES    OF   THREE    MODEL    CYCLONES 

(in  ib) 


TIME 

CYCLONE 

CYCLONE 

CYCLONE 
S 

96 

1011 

1007 

1003 

99 

1008 

1002 

1001 

102 

1006 

1001 

998 

10  5 

1006 

1002 

996 

108 

1005 

1001 

994 

11  1 

1000 

998 

992 

114 

997 

996 

989 

117 

995 

998 

988 

120 

992 

998 

986 

123 

988 

997 

984 

126 

987 

997 

979 

129 

988 

998 

979 

132 

989 

1000 

979 

135 

989 

1001 

978 

138 

987 

1001 

976 

141 

989 

1004 

975 

144 

990 

1005 

975 

147 

988 

1006 

975 

150 

987 

1006 

976 

15  3 

988 

1008 

978 

156 

988 

1008 

980 

159 

985 

1006 

978 

16  2 

985 

1006 

979 

165 

985 

1006 

980 

168 

984 

1006 

981 
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Figure  4. 


Surface  pressure  and  temperature  fields  valid  at 
initial  time  ♦  96  hoars.   Isobars  (solid)  ar«  id 
millibars.   Isotherms  (dashed)  are  in  °c. 
a.  and  b.  as  in  Fig.  3. 
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Figure  5 


10Q0    mb   height    and    10  00   pb.tQ   jiid-layer   H 
thickness    fields   valid   at    initial   time   ♦   96 
hours.      1000    mb    height  contours    {solid)    and 
thickness   contours    (dashed)    are   in   meters. 
a.    and   b.    as    in   Fig.    3. 
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Surface  pressure  and  temperature  fields  at  hour  120  are 
given  in  Fig.  6  with  height  and  thickness  fields  in  Fig.  7. 
During  the  period  from  hour  96  to  hour  123 ,  the  NH  lows  move 
approximately  700  nm  to  the  east  while  moving  only  slightly 
poleward.  Cyclone  N1  decreases  to  992  ab  and  becomes  the 
main  surface  low,  while  cyclone  N2  decreases  to  only  998  mb. 
Approximately  70%  of  the  surface  pressure  decrease  experi- 
enced by  cyclone  N1  occurs  during  the  latter  half  of  this  24 
hour  period.  In  contrast,  more  than  half  of  the  surface 
pressure  decrease  experienced  by  cyclone  N2  occurs  during 
the  first  three  hours  of  this  period.  It  should  be  noted 
that  cyclone  N1  experiences  a  rapid  central  surface  pressure 
decrease  of  17  mb  between  hour  108  and  hour  123.  In  the  SH, 
cyclone  S  decreases  to  986"  mb  at  a  nearly  constant  rate  of 
approximately  two  millibars  per  three  hours. 

The  thermal  ridges  associated  with  cyolones  N1  and  S  are 
more  pronounced  than  at  hour  96;  however,  thermal  advection 
at  the  surface  increases  very  little  during  the  period  from 
hour  96  to  hour  120.  By  the  end  of  this  period,  pockets  of 
warm  air  are  coincident  with  the  low  pressure  centers  in  the 
NH.  However,  cyclone  N2  does*  not  continue  to  develop  and 
thus  will  not  be  discussed  further. 
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Figure  S. 
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Sams    as   Fig.    4    except    for    initial   time   +    120 
hours. 
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Figure  7.      Same   as   Fig.    5    except    for    initial    time   ♦    120 
hours. 
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The  analyses  for  hour  144  are  presented  in  Figs.  8  and 
9.  During  the  period  from  hour  120  to  hour  144,  cyclone  N1 
moves  approximately  700  nm  to  the  east-northeast  (ENE)  while 
cyclone  S  moves  approximately  5  00  nm  to  the  east-southeast 
(ESE).  Cyclone  N1  experiences  little  intensification  as  the 
central  surface  fluctuates  between  992  and  987  mb  with  the 
lowest  value  occurring  at  hours  126  and  138.  The  low  pres- 
sure area,  initially  to  the  northeast  of  oyclone  N1,  inten- 
sifies slightly  and  moves  to  a  location  approximatey  1000  nm 
to  the  northwest  of  cyclone  N1  in  an  area  of  lower-tropos- 
pheric  warm  advection.  During  this  same  period,  the  central 
surface  pressure  of  cyclone  S  decreases  to  975  mb  with  more 
than  half  of  this  decrease  occurring  during  the  first  three 
hours  of  the  period.  At  hour  1 32  (not  shown) ,  a  pocket  of 
warm  air  is  coincident  with  the  surface  center  of  cyclone  S. 
A  significant  feature  of  the  NH  thickness  field  is  the  north 
to  south  orientation  and  the  amplitude  of  the  thermal  ridge. 
In  contrast,  the  SH  thermal  ridge  is  oriented  nearly  west  to 
east   and    is   extremely   elongated. 

Surface  pressure  and  temperature  fields  valid  at  hour 
168  are  presented  in  Fig.  10  with  the  height  and  thickness 
fields  given    in    Fig.    11.      During    the    period    from   hour    144    to 
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Same   as   Fig.    4    except    for   initial    time   ♦   144 
hours. 
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Figure  9.      Same    as   Fig.    5    except    for    initial   time   +   144 
hours. 
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hour  156  hours,  cyclone  N1  remains  relatively  steady-state 
with  a  central  surface  pressure  of  approximately  988  mb 
while  the  central  surface  pressure  of  cyclone  S  rises  5  mb 
to  980  mb.  From  hour  156  to  hour  168,  the  central  surface 
pressure  of  cyclone  N1  decreases  to  984  nb  while  cyclone  S 
remains  nearly  constant  at  about  980  mb.  During  the  final 
24  hcurs  of  this  study,  the  MH  thermal  ridge  shifts  from  a 
north-south  orientation  to  a  northwest-southeast  orienta- 
tion. However,  the  SH  thermal  ridge  maintains  its  west-east 
orientation.  The  1000  mb  height  contours  nearly  parallel 
the  thickness  contours  poleward  of  each  cyclone  so  there  is 
little  warm  advection  in  the  vicinity  of  these  cyclones; 
however,  cold  air  advection  continues  eguatorward  of  each 
cyclone.  During  the      period   from     hour    144      to   hour      162, 

cyclone      N1    moves     approximately      300   nm      to      the   NNE.  It 

appears  that  cyclone  N1  moves  approximately  360  nm  during 
the  last  six  hours  of  this  study;  however,  this  apparent 
rapid  movement  may  actually  be  a  reformation  of  the  low  in 
the  vicinity  of  the  low  pressure  area  located  northwest  of 
cyclone  N1  at  hour  144.  Cyclone  S  also  exhibits  an  apparent 
rapid  movement  during  the  final  24  hours  of  simulation,  mov- 
ing   approximately   720   nm   to     the    ESE.         Approximately    500    nm 
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of  this  movement  occurs  during  the  period  from  hour  150  to 
hour  162  as  the  central  surface  pressure  increases  three 
millibars  to  979  mb.  As  in  the  NH,  it  is  possible  that  this 
apparent  rapid  movement  is  a  reformation  of  cyclone  S.  How- 
ever, unlike  the  NH  case,  there  is  no  secondary  low  or  warm 
pocket  in   the  vicinity   of  the  new   location    of   cyclone    S. 

C.        HID- LEVEL    FEATURES 

In  this  section  mid-tr opospheric  features  are  discussed 
using  mid-layer  4  of  the  model  (approximately  500  mb) . 
Kid- layer  4  height  and  temperature  fields  valid  at  hour  96 
are  depicted  in  Fig.  12  with  the  vorticity  and  wind  fields 
given  in  Fig.  13.  A  weak  short  wave  within  nearly  zonal 
flow  is  apparent  in  each  hemisphere.  The  height  and  temper- 
ature fields  are  nearly  in  phase  so  there  is  little  thermal 
advection  in  either  hemisphere.  Developing  cyclones  N1  and 
S  are  located  under  the  short  wave  ridges  in  the  respective 
hemispheres  with  a  weak  vorticity  maximum  located  300  to  400 
nm  to  the  west  of  each  cyclone.  Mi d-tropospheric  wind  flow 
patterns  indicate  weak  cyclonic  vorticity  advection  (CVA) 
over   each    cyclogenesis   area. 
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Figure   10.      Same  as   Fig.    4    excapt    for    initial  time   ♦    168 
hours. 
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Figure   11.      Same  as   Fig.    5    excapt  for   initial   time   +    168 
hours. 
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Figure    12, 


aid-laver   4    height   and  temperature    fields    valid 
at   initial    time  *    96    hours;      Height   contours 
(solid)    are   in    decameters.      Isotherms    (dashed) 
are    m   °C.      a.    and   b.    as    in   Fig.    3. 
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Figure    13. 


Hid-layer  »  absolutf  .Tortjcity    isec-i)    and   wind 
fields   valid   at   initial   time    *    $6    hours. 
Absolute    vorticity    (x    10-s)    contours   are   solid, 
a-    and   b.    as   m  Fig.    3. 
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Hid- layer  4  fields  valid  at  hour  120  are  presented  in 
Figs.    14      and    15.  The   height      and      therial    fields      remain 

nearly  in  phase  with  weak  cold  advection  immediately  to  the 
west  of  the  center  of  cyclone  N1  and  warm  advection  approxi- 
mately 5  degrees  longitude  further  to  the  west.  During  the 
first  24  hours  of  this  study,  the  amplitude  of  the  NH  short 
wave  doubles  while  that  of  the  SH  short  wave  increases  by 
approximately  70%.  This  amplitude  intensification  is  nearly 
simultaneous  with  the  most  rapid  surface  development  which 
occurs  during  the   period   from   hour  96   to   hour    126. 

The  orientation  of  the  NH  short  wave  trough  is  northwest 
to  southeast  while  that  of  the  SH  short  wave  is  southwest  to 
northeast  as  each  cyclone  maintains  its  position  relative 
its    associated   short    wave.  The  cyclonic:    vorticity  maximum 

in  the  NH  increases  in  magnitude  and  is  located  directly 
over  cyclone  N1  resulting  in  anticyclonic  vorticity  advec- 
tion (AVA)  west  of  the  center  of  cyclone  N1.  The  increase 
in  the  vorticity  maximum  at  mid-layer  4  and  the  most  rapid 
surface  development  also  occur  simultaneously.  In  the  SH, 
the  cyclonic  vorticity  maximum  also  increases  in  magnitude; 
however,  in  contrast  to  the  NH  oase,  the  SH  vorticity  max- 
imum     does      not      move  over      the      storm   center.         Instead   the 
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vorticity  maximum  moves  to  the  east  as  cyclone  5  moves 
southeasterly  to  a  position  poleward  of  the  vorticity  max- 
imum. As  in  the  NH  case,  intensification  of  the  vorticity 
maxiium  occurs  simultaneously  with  the  most  rapid  surface 
development.  Zero  vorticity  advection  occurs  over  cyclone  s 
while   strong   CVA   occurs   to   the   north   and   r.ortheast. 

Mid-tropospheric  fields  valid  at  hour  144  are  given  in 
Figs.  16  and  17.  The  amplitude  of  each  short  wave  continues 
to  increase  and  the  orientation  of  the  trough  axes  becomes 
lore  west-east.  Extremely  sheared  flow  is  evident  in  the  SH 
and,  to  a  lesser  degree,  in  the  NH.  The  short  waves  and 
their  associated  cyclones  maintain  their  positions  relative 
to  each  other  during  the  period  between  hours  120  and  144  as 
the  height  and  thermal  fields  remain  nearly  in  phase. 
Although  maintaining  its  position  over  cyclone  N1,  the  vor- 
ticity maximum  in  the  NH  weakens,  having  reached  a  peak  at 
hour  132  (not  shown) .  As  in  the  NH  case,  the  SH  vorticity 
maximum  reaches  its  peak  at  hour  132  and  naintains  its  posi- 
tion relative  to  cyclone  S.  In  each  case,  the  vorticity 
patterns  are  quite  elongated,  because  of  the  extreme  shear 
present   in   the   wind   field. 
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Figure   14.      same   as  Fig.    12   except   for   initial   time   *■    120 
hours. 
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Pigure   15.      Same   as   Fig.    13  except   for    initial  time  *    120 
hours. 
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Figure   16.      Same   as   Pig.    12  except   for    initial  time  ♦    14U 
hours. 
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Figure   17.      same  as  Fig.    13   except  for   initial   time   *    1UU 
hours. 
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At  hour  168  (Figs.  18  and  19) ,  cyclone  N1  is  positioned 
under  the      inflection   point      ahead  of      the    trough-  The    NH 

short  wave  continues  to  amplify  and  tha  shearing  is  more 
pronounced  as  the  axis  of  the  short  wave  is  nearly  east  to 
west.  Between  hours  144  and  15 8r  the  SH  upper  trough  fills 
leaving  a  split  flow  pattern  with  nearly  zonal  flow  to  the 
south  and  the  weakening  trough  to  the  aorth  of  46S,  the 
latitude  of  cyclone  S.  A  weak  closed  low  appears  approxi- 
mately 700  nm  to  the  northwest,  however  cyclone  S  fails  to 
develop  into  a  deep  cold  core  vortex.  The  vorticity  fields 
continus  to  weaken  and  become  more  elongated  during  the 
final  24  hours  of  this  study.  Broad,  relatively  uniform 
areas  of  cyclonic  vorticity  ara  located  aquatorward  of  each 
cyclone.  These  vorticity  patterns  nearly  parallel  the  wind 
flow  indicating  little  vorticity  advection  in  either 
hemisphere. 

D.        UPPER-LEVEL    FEATURES 

Mid- layer  3  (approximately  3  00  mb)  height  and  tempera- 
ture fields  valid  at  hour  96  are  given  in  Fig.  20  with  the 
wind  fields  given  in  Fig.  21.  Nearly  zonal  flow  is  evident 
with  weak  cyclonic  flow  700  nn  to  the  west  of  each  cyclone 
and    weak    anticyclonic      flow   1000    nm  and   2300    nm      to   the   east 
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of  cyclones  N1  and  S  respectively.  The*  two  jets,  centered 
near  45N  and  35S,  are  nearly  zonal  and  somewhat  stronger 
than  those  in  the  initial  field.  The  NH  and  SH  jet  maxima, 
approximately  45  m  sec-1  and  55  m  sec-1  respectively,  are 
located  900  to  1500  nm  to  the  west  of  the  developing 
cyclones. 

During  the  first  24  hours  of  this  study,  the  flow  imme- 
diately to  the  west  of  each  cyclone  becomes  increasingly 
cyclonic  as  the  amplitude  of  the  trough/ridge  system 
increases.  At  hour  120  (Figs.  22  and  22)  the  cyclones  are 
positioned  slightly  east  and  poleward  of  the  base  of  the 
trough.  The  jet  maximum  in  the  NH  increases  only  slightly 
while  the  SH  jet  maximum  increases  to  over  60  m  sec-1. 
These  jets  are  located  approximately  600  im  west  of  the  base 
of    the  respective   troughs. 

The  upper-level  charts  valid  at  144  hours  are  shown  in 
Pigs.  24  and  25.  In  the  NH,  cyclone  N1  is  located  ahead  of 
the  trough  under  nearly  meridional  flow  with  a  jet  maximum 
of  ever  40  m  sec-1  located  approximately  500  nm  to  the 
southeast  and  a  secondary  jet  maximum,  also  greater  than  40 
m  sec-1,  located  near  the  top  of  the  ridge,  approximately 
600    nm  to   the  north   of  cyclone  N1.      In   the    SH,    there   are   two 
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Figure   18. 


Same   as   Pig.    12   except   for   initial   time   +    168 
hours. 
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Figure   19.      Same  as  Fig.    13   except   for   initial  time  ♦    168 
hours. 
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7 layer   3    height   and  temperature    fields    valid 
at   initial   time  ♦    9b .      Height    contours    (solid) 
are    in   decameters.      Isotherms    (dashed)    are    in 
°C.      a.    and   b.    as   in    Fig.    3. 
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8a4- layer. 3  wind  vectors  and  isotachs  valid  at 
initial  time  ♦  96.   Isotachs  (solid)  in  n  sec-1, 
a.  and  b.  as  in  Fig.  3. 
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Same   as   Fig.    20  except   for   initial  time  +    120 
hours. 
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Figure  23. 
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Same   as   Fig.    21   except  for  initial   time   *■    120 
hours. 
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jet  laxima  (each  greater  than  5  0m  sec-i)  associated  with 
cyclone  S--one  behind  the  trough  and  one  in  the  base  of  the 
trough.  In  each  hemisphere  the  flow  ahead  of  the  trough  is 
basically  meridional  and  the  short  waves  are  continuing  to 
increase    in  amplitude. 

The  mid-layer  3  wind  field  valid  at  hour  156  (not  shown) 
is  characterized  by  several  small  disjointed  jet  maxima 
greater  than  U0  m  sec-1  in  the  NH.  There  appears  to  be  no 
organization  to  these  areas  as  some  ace  ahead  and  some 
behind  the  trough.  By  this  time,  the  SH  trough  has  tilted 
back  on  itself  and  is  oriented  nearly  west  to  east  and  the 
flow  south  of  the  trough  is  now  zonal.  A  developing  cut-off 
low  is  indicated  by  the  extra  me  shear  in  the  SH  height 
field.  Similar  directional  shear  is  prasent  in  the  NH,  alt- 
hough not    to  the   extent   of   that    in  the  SH. 

Figs.  26  and  27  depict  the  mid-layer  3  patterns  at  hour 
168.  In  the  NHr  the  position  of  cyclone  N1  relative  to  the 
trough  remains  unchanged  during  the  final  24  hours  of  this 
study.  However,  an  organized  jat  maximum  of  over  U0  m  sec-1 
to  the  west  is  approaching  tha  base  of  the  trough.  This 
maximum  is  a  result  of  the  cyclic  continuity  assumption 
employed    by   the      model   and    is   the   same   maximum      which    was    in 


59 


a. 


q  h  i  i  1  i  i  i  1  i  i  i  I  i  i  i  1  i  i  i  I  i  i  i  I  i  i  i  I  t  i  i  1  i  i  i  I  i  i  i  1  i  i  i  I  i  i  r 

0  20  40  60  80  100  120 

L0NGITUDE 


Figure  24.      same   as   Fig.    20  except   for   initial  *ime  ♦    144 
hours. 
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Same  as  Fig.    21   except  for   initial   time   «■    144 
hours. 
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the  downstream  ridge  (north  of  the  cyclons)  at  approximately 
hour  144.  The  flow  directly  over  and  north  of  the  cyclone 
N1  remains  basically  meridional  but  shifts  somewhat  to  the 
northwest  producing  stronger  directional  shear  to  the  north- 
west of  cyclone  N1 .  In  the  SH,  a  jet  maximum  greater  than 
60  m  sec-1  is  centered  in  the  base  of  the  trough  750  nm 
equatorward  of  cyclone  S,  which  is  located  under  a  region  of 
very  weak  flow  at  mid-layer  3.  A  weak  closed  low  at  mid- 
layer  3  is  located  approximately  700  nm  northwest  of  cyclone 
S,  almost  directly  over  the  weak  closed  low  at  mid-layer  4. 
The    flow    equatorward   of   each   cyclone   continues   to  be   zonal. 

E.        COHPARISON    TO    OBSERVED    CYCLONES 

A  superficial  analysis  of  model  cyclogenesis  might  give 
the  impression  that  these  storms  are  realistic  in  every 
aspect.  A  closer  examination  reveals  thsm  to  have  a  number 
of    unusual    features.  In    this    section,      the      realistic  and 

anomalous    features    will    be    summarized. 

Although  cyclone  S  is  more  intense,  the  overall 
development  of  both  model  cyclones  resembles  that  of  storms 
observed  in  the  North  Atlantic  by  Petterssen,  et  al.  (1962). 
The  model  cyclongenesis  begins  in  the  lower  troposphere  in 
the    vicinity   of    a   baroclinic     zone,      under   nearly   zonal   flow 
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Figure  26.      Same   as   Pig.    20  except   for   initial  time  ♦    168 
hours. 
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Figure  27. 


Same  as   Fig.    21   except   for   initial   time   +    168 
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aloft,  and  with  no  significant  vorticity  advection.  The 
upper-level  cold  trough  and  the  low-level  cyclone  develop 
simultaneously  while  maintaining  their  position  relative  to 
each  other.  In  addition,  the  primary  contribution  to 
development  is  provided  by  lower -tropospharic  thermal  advec- 
tion (at  least  during  the  early  stages  of  development)  and 
not    upper- level    vorticity   advection. 

During  the  early  stages  of  cyclogenesis,  the  lower- 
tropospheric  features  are  very  realistic.  Near  the  surface, 
the  simulated  cyclones  resemble  the  classical  cyclone  model 
with  the  cold  air  on  the  west  side  of  the  cyclone  and  warm 
air  on  the  east  side.  Cyclones  N1  and  S  are  nearly  symmet- 
rical for  the  first  24  to  48  hours  of  this  study.  However, 
the  subsequent  development  results  in  a  very  elongated 
structure.  although  asymmetrical  cyclones  are  not  uncommon, 
the  extreme  elongation  of  these  cyclones  does  appear  to  be 
atypical . 

The  lo wer-tropospheric  thermal  patterns  are  also  real- 
istic during  the  first  half  of  the  period  studied.  The  warm 
advection  in  the  vicinity  of  the  storm  during  the  early 
stages  of  development  and  the  warm  air  coincident  with  the 
storm   center     during   the   latter      stages   of      simulation    would 
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tend     to      support     cyclogenesis.  However,      the     elongated 

thickness  patterns  and  the  oriantation  of  the  thermal  ridges 
in    each   hemisphere  during   the   later   stages   are   unusual. 

The  Nitta  and  Yamamoto  (1974)  observational  study  sup- 
ports the  findings  of  Petterssen,  et  al.  (1962).  In  addi- 
tion, they  noted  that  oceanic  cyclones  generally  predominate 
in  the  lower  troposphere.  The  development  of  the  two 
cyclones  used  in  this  study  is  generally  limited  to  the 
lower  troposphere  although  cyclone  S  developes  weak  closed 
lows  at  mid-  and  upper-troposphe ric  levels.  These  features, 
however,  are  extremely  weak  and  do  not  appear  until 
hour  168.  It  is  not  surprising  that  cyclone  S  developes 
to  the  upper  troposphere,  however  the  vertical  coupling 
of  this  development  is  somewhat  unusual.  This  feature  will 
be    discussed   in    more   detail  later   in    this   section. 

The  above  characteristics  of  the  nodel  storms  which 
appear  to  be  realistic  are  basically  associated  with  the 
lower  troposphere.  On  the  other  hand,  the  anomalous  fea- 
tures are  generally  associated  with  the  vertical  configur- 
ation of  the  cyclones  and  the  structure  of  the  mid-  and 
upper-tr opospher ic  levels  of  the  model,  mid-layer  4  and 
mid-layer      3   respectively.  The   development      of  the      short 
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waves  in  the  middle  and  upper  troposphere  appears  to  be 
atypical.  In  particular,  the  Drientation  of  these  waves  and 
the  resulting  shear  are  quite  unusual.  as  noted  earlier, 
troughs  in  each  hemisphere  eventually  tilt  back  on 
themselves  creating  extreme  directional  shear. 

The  height  and  thermal  fields  at  mid-layers  4  and  3  are 
nearly  in  phase  throughout  the  7  2  hours  of  this  study.  Fur- 
thermore, a  close  examination  of  the  height  fields  at  these 
two  levels  reveals  them  to  be  nearly  identical  in  shape  and 
vertically  stacked  throughout  the  entire  period  of  this 
study.  This  is  a  reflection  of  a  self-amplification  process 
and  is  likely  related  to  the  lack  of  occlusion  and  shallow 
structure  of  the  mature  cyclones.  In  addition,  the  upper 
troughs  and  the  cyclones  maintain  their  positions  relative 
to  each  other  throughout  the  period  of  this  study.  These 
two  features  are  consistent  with  tha  observations  of 
Petterssen,  et  al.  (1962);  however,  the  distance  between 
cyclone  s  and  the  weak  upper  lows  which  develop  in  the  SH 
after  16  8  hours  of  simulation  is  unusual.  The  closed  upper 
lows  are  vertically  aligned  but  they  are  nearly  700  nm  to 
the  northwest  of  the  surface  center.  By  the  time  the  upper 
lows  develop,   cyclone  S  has  reached  its  peak  intensity  (975 
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mb    between   hours    141    and   1  47)      at. the   surface      and   has   begun 
to      weaken    somewhat      (to      981    mb)  .  One      might    expect      the 

upper-level  lows   to     be    more   intense   and      vertically   aligned 
with  the    surface  low   during  tha    mature   stage   of   the  cyclone. 

The  development  of  the  pattarns  in  the  mid-layer  4  vor- 
ticity  fields  also  appears  to  be  anomalous.  In  particular, 
the  location  of  the  vorticity  maxima  with  respect  to  the 
cyclones  (over  cyclone  N1  and  eguatorward  of  cyclone  S) 
during  the   final   43   hours     of      this   study   is    unusual.  The 

extremely   elongated   regions  of      relatively    uniform   vorticity 
which  parallel   the   wind    flow  also  appear  to   be   unusual. 

In  conclusion,  these  model  storms  are  characterized  by 
both  realistic  and  unusual  features.  The  realistic  features 
are  generally  in  the  lower  troposphere  with  the  atypical 
features  generally  in  the  middle  and  upper  troposphere.  The 
vertical  resolution  of  the  model's  upper  troposphere  could 
possibly  be  a  factor  in  the  development  of  these  unusual 
features.  Although  all  aspects  of  these  storms  are  not  typ- 
ical of  oceanic  cyclones,  they  are  realistic  enough  to  be 
useful  in  determining  the  mass  and  circulation  budgets  of 
numerically   simulated   cyclones.      Quasi-Lagrangian   diagnostic 
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techniques  are  described  in  the  next  chapter.  An  analysis 
of  the  mass  budgets  of  thsse  storms  is  presented  in 
Chapter  V. 
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IV.       20MIrLAGRANGIAN    DIAGNOSTIC    TECHNIQUE 

A.        QUASI-LAGRANGIAN    DIAGNOSTICS 

A  primary  goal  of  this  thesis  is  to  describe  cyclone 
development  quantitatively  through  an  analysis  of  mass  bud- 
gets. The  quasi-Lagrangian  diagnostics  approach  provides  a 
means  of  achieving  this  goal.  Wash  (1973)  demonstrated  the 
validity  of  the  quasi-Lagrangian  diagnostics  approach  in 
budget  studies  of  observed  and  numerically-predicted 
cyclones.  Basic  quasi-Lagrangian  diagnostics  techniques  are 
presented  in  this  section  with  the  mass  budget  equation  pre- 
sented in  the  following  section.  A  description  of  the 
generalized  quasi-Lagrangian  diagnostics  budget  equation  and 
the  budget  volume  coordinate  system  (Wash,  1978)  are 
presented   in   Appendix   B. 

Traditionally,  budget  studies  have  employed  a  stationary 
budget  volume  which  is  positioned  so  that  the  feature  being 
studied  moved  through  or  remained  within  the  budget  volume 
during  the  period  of  the  study  (e.g.  Smith  (1973),  Kung  and 
Baker  (1975),  and  others).  Since  the  horizontal  extent  of  a 
stationary   budget   volume      is   greater    than    that      of   a   cyclone 
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scale  feature,  it  is  not  possible  to  isolate  a  cyclone  from 
the  environment  in  which  it  is  moving.  In  addition,  budget 
studies  employing  stationary  budget  volumes  provide  a  quali- 
tative estimate  of  the  large  scale  forcing  processes 
involved  in  cyclone  development  and  decay.  Therefore,  the 
use  of  a  stationary  budget  volume  is  not  particularly  suited 
to    budget    studies   of   cyclones. 

The  quasi-Lagrangian  technique  is  ideal  for  budget  stud- 
ies cf  cyclones  because  the  budget  volume  moves  with  the 
velocity  and  direction  of  the  cyclone.  A  translating  budget 
volume  affords  a  method  of  distinguishing  between  the 
effects  of  storm  translation  and  actual  cyclone  development. 
Thus  a  translating  budget  volume  is,  in  fact,  an  instrument 
for  examining  the  fluxes  across  the  budget  volume  boundaries 
by  eliminating  the  distortion  that  results  from  the  movement 
of    the  storm. 

Eudget  studies  which  employ  a  storn-f ollowing  budget 
volume  provide  a  view  of  the  interaction  between  the  smaller 
scale  cyclone  and  the  planetary  scale  environment  in  which 
it    is  embedded.  Through    the   transport   relation     and   other 

governing  equations,  the  quasi-Lagrangian  budget  approach 
provides      a      method     of     calculating      the      underlying      scale 
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interaction  processes.  This  approach  emphasizes  the 
transports  of  mass  through  the  boundaries  and  the  changes  in 
the  circulation  within  the  budget  volume.  Quasi-Lagrangian 
diagnostics  incorporates  transports  into  and  out  of  the 
budget  volume,  as  well  as  internal  sources  and  sinks  of 
these   basic   properties. 

The   quasi-Lagrangian   diagnostics   approach      has   been   used 
in      a   variety     of  cyclone      studies.  Table    II,        partially 

extracted  from  Hash  (1978),  is  an  updatad  list  of  several 
rscent  budget  studies  which  applied  quasi-Lagrangian 
diagnostics   to   observed    and/or   simulated  cyclones. 


TABLE   III 

QUASI-LAGRANGIAN    DIAGNOSTICS    BUDGET    STUDIES 
(partially  extracted    from    »ash,    1978) 


PROPERTY  RESEARCHER  (S) 

Available   Potential  Energy  Spaete    1974 

Circulation  Hash   1975 

Hass  Johnson    and   Downey    1976 

Absolute   Angular   Momentum  Johnson    and   Downey    1976 

Kinetic   Energy  Chen   and    Bosart    1977 

Absolute   Angular  Momentum  Wash    1978 

Mass  Roman    1981 

Circulation/  m  -«Q^ 

Angular   Momentum  Conant    1982 
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As  mentioned  earlier,  quasi -Lagrangim  transport  theory 
provides  a  diagnostic  framework  which  is  ideally  suited  to 
the  study  of  cyclone-scale  circulations.  In  addition,  the 
internally  consistent  "data"  base  of  numerically  simulated 
cyclones  limits  computational  arrors  to  those  resulting  from 
interpolation  from  a  rectangular  grid  to  the  circular  grid 
of      the      budget    volume   and   to   round- off  errors.  With      the 

guasi-Lagrangian  diagnostics  approach,  a  cylindrical  budget 
volume  is  placed  around  the  model  storm  and  is  moved  with 
the  storm  center.  For  this  study,  tha  surface  pressure 
minimum  is  used  to  indicate  the  storm  canter.  The  budget 
volume  extends  from  the  earth's  surface  to  the  top  of  the 
atmosphere  with  the  axis  of  the  budget  volume  normal  to  the 
earth's  surface.  The  budget  volume  is  actually  a  portion  of 
a  cone  with  its  vertex  at  tha  earth's  center.  However, 
since  the  tropospheric  height  is  much  less  than  tha  radius 
of  the  earth,  the  budget  volume  radius  is  assumed  to  be 
constant  in  the  vertical.  The  budget  volume  for  this  study 
consists  of  the  seven  model  sigma  lavels  as  defined  in 
Chapter  II.  The  radius  of  the  budget  volume  is  defined  in 
degrees      latitude   relative      to      the     storm    center.  Having 

defined      the      location      and      radius     of      tha      budget   volume. 
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the  budget  of  any  basic  storm  properties  (e.g.  mass, 
circulation,  angular  momentum,  kinetic  enargy,  and  moisture) 
can    te  studied. 

B.        EODGET    EQUATIONS 

For  this  study  the  quasi-La grangian  diagnostics  techni- 
que is  applied  to  mass  budgets  of  two  numerically  simulated 
cyclones.  The  model  experiments  defined  vertical  motion  to 
be  zero  at  the  upper  and  lower  boundaries  of  the  atmosphere. 
Thus  vertical  transports  across  these  boundaries  are 
eliminated   in  the   mass   budget  aquations. 

Table  IV  contains  the  mass  budget  aquation  in  the  sigma 
coordinates  of  the  model.  In  its  most  basic  form,  the  gen- 
eralized budget  equation  consists  of  thrae  terms:  lateral 
transport,    vertical  transport,    and  source    (or   sink). 

Since  there  are  no  sources  or  sinks  of  mass,  the  source 
(sink)         term  of      the   generalized     equation      is   zero.  The 

remaining  terms  in  the  mass  budget  equation  represent  the 
lateral  mass  transport  (convergence/divergence)  and  vertical 
transport  (residual)  respectively.  When  the  limits  of  inte- 
gration are  a  =  0  and  o  -  1  (whare  (7=0),  the  budget  resid- 
ual term  must  be  zero.  If  not,  the  lataral  mass  transport 
term   for   the      various   model  layers   are   adjusted      so   that   the 
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integrated  residual  vanishes.  Lateral  mass  transport  and 
mass  tendency  are  computed  within  the  mass  budget  and  resi- 
dual imbalances  are  then  computed  from  the  time  averages  of 
these  quantities.  Vertical  integration  of  the  layer  resid- 
uals produces  a  vertical  mass  transport  which  is  then 
expressed   as   a   vertical   velocity    field. 


TABLE    IV 
QUASI-LAGRANGIAN    HASS    BUDGET    EQUATION    IN    SIGMA    COORDINATES 

(after  Wash,    1978) 

(See  Appensix   B   for   budget  notation 
and   generalized   budget   equation.) 


For   the  mass  budget   f  ■   1  and  gr        ■     0« 


Substitution   of  these  values   into  the 
generalized   budget   aquation  results   in 
the   following   mass    budget   equation: 


a.T      2ir 


ff  -  J         J   -y  (U-W)6  r  sin6da(-da)f8B 
QB         0 


Q-T       ?B  2TT 


+ 
"a 


2l    |£L  r2  sinBdadB(-da) 
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V.       BUDGET    AN  ALT  SIS 

A.        INTSODOCTION 

The  analysis  of  the  guasi-Lagrangian  mass  budget  begins 
with  a  presentation  of  the  thermal  structure  of  the  cyclones 
in  terms  of  potential  temperature  (Section  B)  .  This  is 
followed  by  a  discussion  of  the  pressure  (mass)  tendency 
(Section   C)  and  the       total       lateral      mass     transport 

(Section    D).  Vertical      mass      transports     are     discussed 

relative  to  the  model  vertical  velocities  (Section  E) .  In 
an  effort  to  determine  the  accuracy  of  the  vertical  velocit- 
ies computed  in  the  mass  budget,  budget  vertical  velocities 
are    compared   to    the    model   vertical  velocities    (Section   P) . 

Host  figures  presented  in  this  chapter  are  layer/time 
pericd  cross-sections.  Budget  calculation  values  at  levels 
and  "data"  times  were  average!  to  produce  layer  and  time 
period  midpoint  values  respectively  for  the  time  sections. 
Exceptions  are  indicated  where  appropriate.  Inward  and 
upward  transports  are  represented  by  positive  values  in  all 
these   figures. 
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B.        THERMAL    STRUCTURE 

Area-averaged  mid-layer  and  "data"  time  potential  temp- 
eratures are  presented  in  Figs.  28  and  29.  The  vertical 
structures  are  similar  in  the  two  cyclones.  The  patterns 
for  radii  six  and  nine  are  nearly  identical.  In  the  upper 
troposphere,  the  potential  temperature  increases  slightly, 
while  decreases  are  seen  in  the  middle  and  lower  troposphere 
during  the  last  half  of  the  study.  During  the  period  of 
most  intense  cyclogenesis  (hours  96  to  125)  only  small  fluc- 
tuations are  seen  in  the  area-averaged  potential  tempera- 
ture. The  cooling,  which  occurs  during  a  quasi-steady  state 
pericd  for  each  cyclone,  is  more  pronounced  in  cyclone  N1. 
This  cooling  within  each  budget  volume  is  a  consequence  of 
the  development  and  poleward  movement  of  the  cyclones.  Tn 
addition,  although  the  troposphere  of  cyclone  N1  cools  more 
than  cyclone  S,  its  average  potential  temperature  is  only 
slightly  lower  at  the  final  "data"  time,  since  cyclone  N1 
was  initially  warmer  in  the  low  troposphere.  These 
potential  temperature  cross-sections  show  that  cyclone  N1 
initially  had  weaker  vertical  stability  than  cyclone  S. 
However,  by  the  end  of  the  study,  the  two  cyclones  have 
similar   vertical   stability   characteristics. 
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The  cooling  in  cyclone  N1  gives  the  appearance  of  a  cold 
core  structure.  This  however  is  not  the  case  since  the  cen- 
ter of  each  cyclone  is  actually  located  in  the  warmest  air 
within  each  budget  volume  (see  Fig,  11).  The  cooling  is 
likely  a  result  of  cold  advection  into  the  eguatorward  half 
of  each  budget  volume.  The  low- tropospheric  thermal  pattern 
can  be  seen  in  the  1000  mb  to  mid-layer  %  thickness  pattern 
fcr  hour  168  (see  Fig.  11).  The  thickness  pattern  shows  a 
warm  tongue  extending  from  the  east  through  the  center  of 
cyclone  N1 ,  while  a  similar  warm  tongue  extends  east  to  west 
through  the  radius  nine  budget  volume  of  cyclone  S.  The 
thickness  pattern  also  reveals  that  the  cold  air  occupies 
more  of  the  budget  volume  of  cyclone  N1  than  cyclone  S,  thus 
resulting  in  the  slightly  lower  average  temeperatures  for 
cyclcne  II. 

Although  the  model  atmosphere  has  no  specified  tropo- 
pause, Figs.  28  and  29  clearly  show  the  existence  of  a 
strong  tropopause  and  a  stratosphere  during  the  72  hours  of 
this  simulation  study.  Level  2  (approximately  200  mb) 
effectively  functions  similar  to  a  tropopause  layer  with 
slight  warming  and  strong  vertical  stability  above  200  mb. 
The  potential   temperature  evolutions   are  characterized   by 
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slight      fluctuations      below      level  2.  These      fluctuations 

provide  a  hint  of  a  cyclic  evolution  pattern  during  the 
development  of  these  storms,  although  the  period  of  the 
cycle  is   not   evident. 

C.        HASS    TENDENCY 

There  are  no  sources  or  sinks  of  mass  within  the  atmos- 
phere, however  it  is  possible  for  a  limited  volume  of  atmos- 
phere to  gain  or  lose  mass.  The  time  rate  of  change  of  mass 
within  a  storm- following  budget  volume  provides  a  quantita- 
tive measure  of  cyclogenesis.  Within  the  mass  budget,  the 
total  mass  within  each  specified  budget  volume  is  computed 
for  each  "data"  time.  From  these  values,  the  mass  tendency 
between  "data"  times  is  determined,  and  subsequently 
converted   to   a   pressure   tendency. 

The  mass  tendency  at  3-hour  intervals  for  both  cyclones 
is    presented    in    Fig.    30.  A      cyclic      pattern      is      clearly 

indicated.  These   time   profiles   also      indicate        that      the 

development  of  the  model  cyclones  occurs  in  two  phases--an 
initial  explosive  cyclogenesis  period  (hours  96  to  126)  and 
a  quasi-steady  state  period  (hoars  126  to  168).  A  12-hour 
oscillation  is  clearly  evident  throughout  these  two  phases. 
During  the   explosive      cyclogenesis  phase,      the  cycle      in  the 
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radius  six  budget  consists  of  a  maximum  3-hour  pressure  fall 
followed  by  a  9-hour  period  of  relatively  smaller  surface 
pressure  change.  The  cyclic  pattern  during  the  quasi-steady 
state  phase  is  characterized  by  alternating  periods  of 
deepening  and  filling.  The  12-hour  cycle  can  be  seen  in  the 
pressure  tendency  patterns  for  both  the  inner  and  outer 
radii.  The  pressure  tendency  trace  for  each  budget  volume 
shows  that  thecyclone's  area-averaged  surface  pressure  drops 
in  discrete  intervals  like  a  stsp  function  rather  than  show- 
ing continuous  deepening.  These  oscillations  could  be 
removed  by  using  an  highpass  filter.  Filtering  would  also 
serve  to  emphasize  the  deepening  during  the  first  2U  hours 
and    the   filling    which   occurs   later. 

The  pressure  tendency  oscillations  are  not  as  pronounced 
in  the  graphs  of  the  central  surface  pressures  of  these 
cyclones  (Fig.  31).  This  is  especially  true  during  the 
explosive  cyclogenesis  phase,  when  a  generally  steady 
decrease  is  seen  in  the  surface  pressure.  However,  the  rate 
of  change  nearly  corresponds  to  the  pressure  tendency  fluc- 
tuations in  Fig.  30.  During  the.  second  phase  of  cyclone 
development,  oscillations  are  seen  in  the  central  surface 
pressure,  however  they  are  not  as  large  as  those  the 
pressure   tendency. 
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Area-averaged   pressure  tendency    for   cyclones 
N1   and    S-       (mb/3   hours)      a.    Cyclone    Nl. 
b.    Cyclone   S. 
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The  cause  of  these  mass  fluctuations  is  not  known.  One 
possibility  is  that  the  cycle  is  the  refaction  of  an  exter- 
nal gravity  wave  propagating  within  the  model.  An  external 
wave  would  have  a  similar  effect  on  the  surface  pressures  of 
neighboring  anticyclones.  For  this  reason  a  study  of  the 
central  surface  pressures  of  model-generated  anticyclones 
was    completed. 

This  preliminary  study  of  five  highs  (three  in  the  NH 
and  two  in  the  SH)  for  periods  of  48  to  72  hours  (the  length 
of  the  period  was  dependent  on  whether  the  anticyclone  was 
an  analyzed  feature  at  a  particular  "data"  time)  strongly 
indicates  a  cyclic  pattern  of  surface  pressure  fluctuations 
(approximately  two  millibars)  in  the  moiel  anticyclones  as 
well.  The  pressure  oscillations  of  the  cyclones  and  anticy- 
clones  suggest   a   propagating   external   wave    within   the    model. 

The  pressure  tendency  traces  show  larger  fluctuations 
for  radius  nine  than  radius  six.  This  indicates  that  the 
pressure  falls/rises  for  the  area  between  radii  six  and  nine 
are  larger  than  those  for  the  area  within  the  radius  six 
budget  circle.  An  examination  of  other  budget  radii  showed 
that  the  magnitude  of  the  pressure  teniency  fluctuations 
increases      for      larger        radii.  This        also      suggests      a 
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large-scale  (perhaps  hemispheric  or  global)  sloshing. 
Because  the  scale  of  this  feature  is  luch  larger,  and  it 
does  not  interact  with  the  synoptic-scale  features,  the 
slosh     is      independent      of      cyclone     development.  A      more 

detailed  investigation  of  the  model  surface  pressure  field 
is  required  to  isolate  this  phenomena  and  determine  if  it 
adversely    affects   the   numerical    simulations. 

D.        LATERAL    HASS   TRANSPORT 

One  goal  of  this  thesis  is  to  isolate  the  lateral  and 
vertical  branches  of  the  mass  circulation  of  the  model 
cyclones.  The  lateral  mass  transport  which  is  the  mass  flux 
through   a      specified   budget      volume   boundary.  The   lateral 

mass  transport  time  sections  (Figs.  32  and  33)  graphically 
describe  a  classical  cyclone  structure  of  low-level  converg- 
ence (surface  to  approximately  level  4:  surface  to  600  mb) 
and  upper-level  divergence  (approximately  level  4  to 
lid-layer  2:  600  to  150  mb).  A  similar  pattern  is  noted  at 
radius  nine,  however  the  level  of  nondivergence  (LND)  is 
slightly  higher  (approximately  nid-layer  4:  500  mb)  .  The 
lateral  mass  transport  time  sections  indicate  that  the  mass 
circulation  does  not  extend  into  the  top  two  model  layers. 
This  further  substantiates  the  existance  of  tropopause  and 
stratosphere   behavior  within  the    model. 
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The  maximum  inflow  at  the  inner  and  outer  radii  of  each 
cyclone  occurs  in  the  lowest  model  layer  (centered  at 
approximately  900  mb) ,  while  the  maximum  outflow  is  centered 
near  mid-layer  3  (approximately  300  mb)  at  the  tropopause. 
In  addition,  the  mass  circulation  at  radius  nine  of  each 
cyclone  is  more  intense  during  the  quasi-steady  state  phase 
than  during  the  rapid  cyclogenesis  phase.  This  is  a  further 
indication  that  the  the  second  phase  of  development  is  actu- 
ally an  outward  growth  period  more  than  a  deepening  period 
for    these   cyclones. 

A  cyclic  pattern  is  evident  at  the  levels  of  maximum 
inflow  and  outflow  at  both  the  inner  and  outer  radii  of  each 
cyclone.  Looking  first  at  cyclone  N1  (Fig.  32)  ,  the  maximum 
outflow  at  radius  six  occcurs  at  hour  117  with  secondary 
maxima  centered  at  hours  140,  150,  and  168  while  nearly 
equal  inflow  maxima  are  seen  at  approximately  these  times. 
The  maxima  at  hour  117  occur  luring  the  period  of  intense 
surface  development  (as  indicated  by  changes  in  the  central 
surface  pressure,  Fig.  31)  ,  while  the  maxima  at  later  times 
occur  luring  the  quasi-steady  state  phase  of  cyclone  N1. 
The  lateral  mass  transports  at  radius  nine  of  cyclone  HI 
are        also        characterized      by      oscillations      in      the      basic 
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inflow-outflow     pattern.  Tha      cyclic      variation     in      the 

lateral   mass   transport     are     caapled   with     the     oscillations 
in    the  pressure   field. 


MMt* 


Cycton*  N1 

Cyeioti*  S 


Figure  3  1.   Central  surface  pressures  (mb)  of  cyclones  N1 
and  5. 


A  dramatic  increase  at  6°  and  a  decrease  at  9°  in  the 
height  of  the  level  of  LND  is  also  seen  during  tha  quasi- 
steady  state  phase  of  cyclone  N1  .  The  LND  height  change  for 
radius  nine  (approximately  hour  135)  is  related  to  the 
intensification  of  the  mass  circulatian  during  tha  second 
stage  of  development  and  the  grDwth  of  tha  cyclone  to  the  9° 
budget  volume  boundary. 

The  lateral  transport  tima  sections  for  cyclone  S  also 
exhibit  the  classical  mass  circulation   pattern.    However, 
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the  cyclone  S  lateral  transport  pattern  is  noticeably  dif- 
ferent from  that  of  cyclone  N1.  Specifically,  the  cyclone  S 
lateral  mass  transports  at  radius  six  are  more  intense  than 
those  of  cyclone  VI,  This  is  a  reflection  of  the  greater 
intensity  of  the  spring  hemisphere  cyclone.  In  addition, 
the  height  of  the  LND  at  radius  six  of  cyclone  S  gradually 
lowers  from  about  500  mb  to  near  700  mb,  and  is  generally 
lower  than  its  NH  counterpart  (which  remains  near  600  mb 
throughout  this  study)  after  hour  13  2. 

The  radius  nine  time  sections  for  these  two  cyclones 
also  display  similarities  and  differences.  As  in 
cyclone  VI,  there  is  a  simultaneous  decrease  in  the  height 
of  the  LND  and  the  intensification  of  the  lateral  mass 
transports.  However,  this  occurs  earlier  in  cyclone  S  (hour 
120  for  cyclone  S  and  hour  135  for  cyclone  N1) ,  and  the 
decrease  in  the  LND  height  is  aore  dramatic  (from  500  mb  to 
700  mb  for  cyclone  S  and  about  475  nb  to  600  mb  for 
cyclone  N1).  The  sharp  decrease  in  the  height  of  the  radius 
nine  LND  in  cyclone  S  occurs  simultaneously  with  the  outward 
growth  of  the  cyclone  to  the  outer  budget  volume  boundary. 
As  with  radius  six,  the  radius  nine  lateral  transports  are 
more  intense  in  cyclone  S. 
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Of  special  note  is  the  abrupt  weakening  of  the  cyclone  S 
lass  circulation  during  the  final  21  hours  of  this  study. 
This  weakening  is  most  noticeable  and  occurs  earlier  at 
radius  six.  The  decrease  in  the  mass  transport  intensity 
can  te  seen  in  both  the  inflow  and  outflow  branches,  however 
it  is  more  pronounced  in  the  outflow  branch.  This  weakening 
of  the  lateral  mass  transports  is  indicative  of  the  decay  of 
cyclone  S  and  is  reflected  in  a  corresponding  increase  in 
the  central  surface  pressure  of  cyclone  S.  In  contrast,  the 
lateral  mass  transports  at  the  cyclone  N1  budget  radii  con- 
tinues to  be  strong  at  the  last  "data"  time  when  cyclone  N1 
reaches  its  lowest  central  surface  pressure  (984  mb) . 

E.   VERTICAL  MASS  TRANSPORT 

The  vertical  branch  is  an  integral  part  of  a  cyclone(s 
mass  circulation.  Since  the  vertical  mass  transports  are 
the  area-averaged  vertical  velocity  strictures  within  the 
cyclone  (in  both  isobaric  and  sigma  coordinates) ,  time  sec- 
tions of  the  model  vertical  velocities  (Figs.  34  and  35)  are 
used  for  the  following  discussion  of  the  vertical  trans- 
ports. This  discussion  points  out  similarities  and  differ- 
ences in  the  patterns  of  these  two  storms.  In  general, 
upward  motion   is  seen  below  level  2   (approximately  200  mb) 
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at  the  inner  and  outer  radii.  Alternating  periods  of  weak 
upward  and  downward  motion  are  seen  above  level  2  in  the 
model  stratosphere.  The  upward  motion  is  quite  intense  at 
radios  six  of  each  cyclone.  The  level  of  maximum  upward 
vertical  motion  coincides  with  the  LND  indicated  in  the  lat- 
eral mass  transport  time  sections  (Figs.  32  and  33). 
periods  of  maximum  upward  vertical  motion  coincied  (as  they 
must)  with  the  periods  of  most  intense  Lateral  mass  trans- 
port, transport.  The  upward  vertical  motion  at  radius  nine 
is  weaker  than  at  radius  six.  However,  the  intense  upward 
vertical  motion  occurs  only  within  the  troposphere. 

Considering  first  the  six  degree  budgat  volume,  two  dis- 
tinct periods  of  maximum  upward  vertical  motion  are  evident 
in  each  cyclone.  For  cyclone  N1 ,  the  first  period  (hour  107 
to  123)  occurs  during  the  explosive  cyclogenesis  phase, 
while  the  second  period  (after  hour  147)  occurs  during  the 
quasi-steady  state  phase.  The  subsequent  radius  six  maximum 
occurs  simultaneously  with  a  similar  maximum  at  radius  nine, 
and  results  from  the  outward  growth  of  the  cyclone  to  the 
outer  budget  volume  boundary. 

Cyclone  S  also  experiences  two  periods  of  intense  upward 
vertical  motion   (centered  at   hours  117   and  138).     As  in 
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Fiqure  32.   Cyclone  N1  lateral  mass  transport  tjme  section, 
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cyclone  N1 ,  the  first  maximum  occurs  during  the  explosive 
cyclogenesis  phase  while  the  second  maxiaum  and  its  radius 
nine  counterpart  occur  during  the  period  of  outward  lateral 
growth.  Near  the  beginning  of  the  integration,  a  third 
period  Df  intense  upward  vertical  motion  is  evident  for  each 
budget  radius.  However,  its  significance  is  unclear  due  to 
the  lack  of  the  preceeding  model  data.  In  addition,  the 
decay  of  the  cyclone  is  indicated  by  the  rapid  decrease  in 
upward  vertical  motion  after  hour  159. 

An  interesting  feature  of  the  cyclone  S  radius  nine  time 
section  is  the  downward  vertical  motion  indicated  in  the 
lower  two  layers  between  hours  101  and  111.  This  requires  a 
downward  motion  within  the  5°-9°  area  which  more  than  off- 
sets the  upward  motion  within  radius  six.  This  indicates 
strong  subsidence  around  the  cyclone  compensating  for  the 
intense  upward  vertical  motion  associated  with  the  explosive 
cyclogenesis  phase  when  cyclone  S  is  still  within  the  6° 
budget  volume.  This  is  a  possible  reflection  of  the  cyclic 
evolution  of  the  model  cyclones. 

Time  traces  of  area- averaged  vertical  velocities  for 
level  4  of  cyclones  N1  and  S  are  presented  in  Fig.  36.  The 
two  major   periods  of  strong  upward   vertical  motion   which 
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'igure   34.      Cyclone   N1    model   vart^cal   valpcity   time  section. 
(-//bar   sec-1)       a.    ana  5.    as    in    Fig.    28. 


94 


a. 


'  I  /i  1  I  I  i  i  i ■  ■  ■  i  t  i  I 

»    IOS   111    U7   123   129   135   141   147   153   159   165 

TPC  PERI80S 


b. 


Figure  35.      Cyclone   S   model  vertical   velocity   time  section, 
(-//  bar   sec1)      a.    and  b.    as   in    Fig-    28. 
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correspond  to  the  two-phase  cyclone  development  are  clearly 
depicted.  The  outward  growth  phase  of  each  cyclone  is  viv- 
idly displayed  in  these  traces.  A  dramatic  increase  in  the 
radius  nine  area-averaged  vertical  velocities  is  seen  as  the 
lateral  growth  of  these  cyclones  excseds  the  9°  budget 
volume  boundary. 

The  radius  six  trace  for  cyclone  N  is  characterized  by 
two  distinct  maxima  at  hours  115  and  168.  The  radius  nine 
trace  essentially  has  a  singla  maximum  which  occurs  during 
the    second   phase   of   cyclone  development.  The   extreme   max- 

imum at  hour  115  is  approximately  in  the  middle  of  the  rapid 
cyclogenesis  phase  when  the  central  surface  pressure 
decreases    17    mb    between   hours    13  8   and    123    (see   Fig.    3  1). 

A  surprising  feature  of  the  vertical  velocities  occurs 
between  hours  14  1  and  150  when  the  magnitude  within  radius 
nine  actually  exceeds  the  corresponding  value  for  radius 
six.  This  implies  that  the  vartical  velocities  within  the 
go-go  ragion  exceed  those  within  the  radius  six  budget 
circle. 

Although  the  cyclic  pattern  seen  in  cyclone  S  is  similar 
to  that  of  cyclone  N,  there  are  significant  differences.  As 
in      cyclone        N,        two      maxima        are      seen      at        radius      six 
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(hours  118  and  136),  however  the  time  interval  between  these 
two  cccurrances  is  much  shorter  than  in  cyclone  N.  The 
first  maximum  occurs  during  the  first  phase  of  cyclone 
development  when  the  central  surface  pressure  decreases 
22  mb  between  hours  99  and  126.  The  second  maximum,  which 
occurs  during  the  quasi-steady  state  phase  of  cyclone  S, 
coincides  with  the  single  maximum  at  radius  nine.  The  maxi- 
mum at  radius  nine  of  cyclone  5  is  larger  than  its  cyclone 
N  counterpart.  The  decay  of  cyclone  S  is  clearly  indicated 
by  the  sharp  decrease  in  the  level  4  vertical  velocities  at 
both  six  and  nine  degrees. 

F.   COMPARISON  OF  MODEL  AND  BUDGET  VERTICAL  VELOCITIES 

One  objective  of  this  thesis  is  a  comparison  of  the 
model  vertical  velocities  with  those  computed  in  the  mass 
budget.  The  vertical  divergence  of  the  vertical  mass  trans- 
port (budget  residual)  and  the  mass  tendency  are  computed 
within  the  mass  budget  (see  Tables  IV  and  V).  A  vertical 
integration  of  these  values  (from  an  upper  boundary  value  of 
-  0  at  =0)  describes  the  vertical  mass  transport,  which 
is  then  converted  to  an  area- averaged  vertical  velocity. 
Thus  the  budget  vertical  velocities  are  essentially  kine- 
matic omega  values  for  the  specified  budget  volume. 
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Figure   36.      Model   vertical   velocity  for   level   4. 

(-  fl  bar  sec-1)  a.     ana  b.    as    in    Fig-    30, 
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This  approach  has  been  used  in  atmospheric  studies  to 
infer  vertical  velocities  and  vertical  transports  of  mass 
and  other  properties.  In  this  study,  vertical  velocities 
are  available  as  part  of  the  model  "data"  base.  This 
affords  the  unique  opportunity  to  test  the  accuracy  of  the 
budget-derived  vertical  velocities. 

The  basic  patterns  in  the  lodel  vertical  velocity  time 
sections  (Pigs.  34  and  35)  are  contained  in  the  budget  ver- 
tical velocity  rime  sections  (Figs.  37  and  38) .  The  budget 
vertical  velocities  for  radii  six  and  nine  of  each  cyclone 
are  most  intense  below  the  tropopause  (approximately  mid- 
layer  3)  .  Very  weak  vertical  motion  is  seen  in  the  model 
stratosphere.  At  radius  six  of  each  cyclone,  two  maxima  are 
seen  in  the  budget  vertical  velocities  time  sections,  with 
the  earlier  maximum  occurring  during  the  first  phase  of 
cyclone  development.  As  in  the  model  vertical  motion  pat- 
tern, the  second  maximum  in  the  budget  vertical  velocity 
time  section  occurs  in  conjunction  with  a  radius  nine  maxi- 
mum. This  again  indicates  the  lateral  growth  of  the 
cyclones  beyond  the  9°  budget  volume  boundary.  For  each 
budget  radius  of  cyclone  S  there  is  a  significant  decrease 
in   the  budget   vertical  velocities   after   hour   159   when 
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cyclone  S  begins  to  fill.  This  feature  is  also  seen  in  the 
model  vertical  velocities. 

Although  the  patterns  in  the  model  and  budget  vertical 
velocity  time  sections  are  similar,  significant  differences 
were  noted  in  the  magnitudes  of  the  budget  vertical  veloci- 
ties. Timesections  of  the  differences  between  these  two 
values  (budget  value  minus  modal  value)  are  presented  in 
Figs.  39  and  UO.  With  the  exception  of  the  stratosphere, 
the  budget  vertical  velocities  were  larger  than  those  of  the 
model.  It  appears  that  the  largest  departures  occur  where 
the  vertical  velocities  (both  model  and  budget)  are  largest. 

Since  the  model  has  an  internally  consistent  "data"  base 
there  are  two  possible  sources  for  the  errors  in  the  budget- 
derived  vertical  velocities.  Errors  can  result  from  the 
interpolation  of  the  winds  from  the  rectangular  grid  of  the 
model  to  the  circular  grid  of  the  mass  budget.  Other  errors 
arise  from  the  vertical  distribution  of  the  imbalances 
between  the  integrated  divergence  and  the  mass  tendency, 
which  must  balance  since  there  are  no  sources  or  sinks  of 
mass.  The  comparison  of  the  budget  and  model  vertical 
velocities  indicates  that  the  lateral  mass  transports  are 
too  large  after  the   interpolation  and  balancing  proceedures 
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are  completed  within  the  mass  budget.  Si\ren  the  consistent 
"data"  of  the  model,  these  errors  are  generally  larger  than 
one  night  anticipate.  A  detailed  examination  of  the  inter- 
polation scheme  and  the  vertical  distribution  of  the  budget 
residuals  should  be  conducted  in  an  effort  to  isolate  the 
source  of   these    errors. 
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Figure  37.   Mass  budget  vertical  velocity  time  section  for 
cyclone  s1.   (-/ibar  sec-1) 
a.  and  b.  as  in  Fig.  28. 
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Figure  38.      Mass    budget   vertical    velocity   time   section    for 
cyclone   5.       (7  H  bar    sec-1) 
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a.    and   b.    as  'in  Fig.    28. 
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Figure  39.   Budaet  vertical  velocity  error  time. section  for 
cyclone  N1.   Model  vertical  velocities  minus 
budget  vertical  velocities.   (-/ibar  sec-1) 
a.  and  b.  as  in  Fiq.  28. 
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Figure    40. 


Budget  vertical  velocity  error  time  section  for 
cvclone  S.   Model  vertical  velocities  minus 
tudget  vertical  velocities.   (-/ibar  sec-1) 
a.  and  b.  as  in  Fig.  28. 


105 


VI.       CONCLUSIONS 

In  this  thesis  two  numerically  simulated  oceanic  extra- 
tropical  cyclones  were  studied.  Initially  a  synoptic  study 
of  these  cyclones,  including  a  comparison  with  observed 
cyclones,  was  completed.  These  cyclones  were  found  to  be 
characterized  by  both  realistic  and  unusual  synoptic  fea- 
tures. In  a  general  sense,  the  realistic  features  were 
associated  with  the  lower  troposphere  while  the  unusual  fea- 
tures were  associated  with  the  middle  and  upper  troposphere 
and    the    vertical  coupling   of   the    model   layers. 

These  cyclones  also  were  studied  through  the  application 
of  quasi-Lagrangian  diagnostics.  This  diagnostic  study  of 
the  mass  circulation  of  model  oyclones  showed  a  concentra- 
tion of  inward  lateral  transport  (convergence)  in  the  lowest 
model  layer  and  outward  lateral  transport  (divergence)  in  a 
layer  centered  near  300  mb.  Although  the  model  has  no  spec- 
ified tropopause,  the  basic  patterns  of  the  latsral  and 
vertical  mass  transports  and  the  thermal  structure  revealed 
the    existence  of  a   tropopause   near  mid-layer   3   of  the    model. 

Traces  of  the  mass  tendency  within  6°  and  9°  budget 
volumes      revealed     a      12-hour   oyclic   pattsrn      imbedded    in      a 
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two- phase  cyclogenesis  pattern.  The  first  phase  of  cyclone 
development  was  a  period  of  ax  plosive  dsepening  while  the 
second  phase  was  a  quasi-steady  state  period.  The  12-hour 
cycle  was  most  easily  seen  in  the  pressure  tendency  trace. 
During  the  first  phase  of  cyclone  development,  the  cycle 
consisted  of  a  maximum  3-hour  pressure  fall  followed  by  a 
9- hour  period  of  little  surface  pressure  change.  During  the 
second  phase,  the  cycle  was  characterized  by  alternating 
periods  of  pressure  falls  and  rises.  It  was  also  noted  that 
the  fluctuations  in  the  pressure  tendency  increased  for  suc- 
cessively larger  budget  radii.  The  two-phase  cyclone 
development  was  seen  in  the  traces  of  the  central  surface 
pressures.  However,  the  12-hour  cyclic  pattern  was  not  eas- 
ily detected  in  these  traces.  This  12-hour  cycle  produces  a 
step  function  decrease  in  the  area-averaged  surface  pressure 
of  these  cyclones  rather  than  a  steady  deepening.  This  step 
function  evolution  is  not  clearly  indicated  by  the  central 
surface  pressure  traces  (see  Fig.  31).  A  possible  cause  of 
this  cycle  is  an  external  wave  propagating  within  the  model. 
However,  further  studies  are  needed  to  isolate  the  source  of 
these  oscillations. 
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The  lateral  mass  transport  and  the  vertical  mass  trans- 
port (converted  to  a  vertical  velocity)  clearly  described 
the  mass  circulation  of  these  cyclones.  Time  sections  of 
these  terms  also  indicated  two-phase  cyclone  develpment. 
Two  maxima  were  seen  at  radius  six  of  each  cyclone — one 
during  each  phase  of  development.  The  one  maximum  seen  at 
radius  nine  occurred  in  conjunction  with  the  second  maximum 
at  radius  six  and  indicated  the  outward  growth  of  the 
cyclones  beyond  the  outer  budget  volume  boundary.  Maxima  in 
the  lateral  mass  transport  and  vertical  velocity  time  sec- 
tions generally  occurred  simultaneously.  In  addition,  the 
onset  of  the  decay  of  cyclone  S  was  indicated  by  the  rapid 
decrease  in  the  intensity  of  the  lateral  mass  tendancy  and 
vertical    motion    near   the   end   of    the  study. 

Area-averaged  vertical  velocities  were  computed  within 
the  mass  budget  and  were  compared  to  area-averaged  model 
vertical    velocities.  Although    tha    time    sections     of    these 

terms  have  similar  basic  pattarns,  significant  differences 
were  found  in  their  magnitudes.  It  appeared  that  the  depar- 
tures were  largest  where  the  vertical  velocity  (both  model 
and  budget)  were  the  largest.  Tha  two  possible  sources  of 
these   errors     are    interpolation    of      the   wind   from      the    model 
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grid  to  the  budget  volume  grid,  and  imbalances  in  the 
vertical  distribution  of  the  mass  budget  residuals.  Further 
studies  are  required  to  isolate  the  relative  magnitudes  of 
these  two   errors. 
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711.       RECOMMENDATIONS    FOR    FURTHER    STUDIES 

The  results  of  the  synoptic  study  and  the  mass  budget  analy- 
sis suggest  that  further  diagnostic  studies  of  the  basic 
properties  of  these  cyclones  'are  needed  to  fully  understand 
their  dynamic  structure.  In  addition,  further  studies  are 
needed  to  determine  the  source  of  the  errors  in  the  budget 
derived  vertical  velocities.  With  this  in  mind,  the  follow- 
ing   recommendations   are   made: 

•  Examine  the  budget  technique  for  computing  area- 
averaged  vertical  velocities  and  the  interpolation 
scheme  employed  within  ths  mass  budget  to  isolate  the 
sources  of  the  errors  in  the  budged  derived  vertical 
velocities. 

•  Perform  budget  studies  on  the  model  cyclones  during 
earlier  and  later  time  periods  than  covered  in  this 
study. 

•  Perform  budget  studies  on  model  generated  anticyclones 
to  determine  if  they  are  characterized  by  the  *  2-hour 
cyclic  pattern   of   the    modal    cyclones. 

•  Study  the  cyclic  oscillations  by  differencing 
successive    pressure   fields. 

•  Perform  budget  studies  on  cyclones  and  anticyclones 
generated  in  the  Sandgathe  (1981)  adiabatic 
experiments. 

•  Perform  budget  studies  of  other  oroperties  on  these 
model   cyclones   and    on    model    anticyclones. 
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APPENDIX    A 
COMMONLY    OSED    ACRONYMS    AND    ABBREVIATIONS 

N03APS        Navy  Operational  Atmospheric   Prediction    Model 
UCLA        University   of  California,    Los    Angeles 


NNW        Naval   Environmental   Prediction    Research 
Facility 


ENE  east-northeast 

ESE  east-southeast 

NNE  north- northeast 

SSI  south-southwest 

LND  level   of   nondivergence 

AVA  anticyclonic   vorticity   advection 

CVA  cyclonic   vorticity    advection 

PBL  planetary   boundary    layer 

see-1         per   second 

mb  millibar 

°C  degree   Celcius 

°K  degree    Kelvin 

1 11 


NH  northern   hemisphere^ 

SH  southern   hemisphere 

no  nautical   miles 

H  mass 

N  north 

S  south 

on  meters 
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APPENDIX    B 

QOASI-LAGRANGIAN   30DGET    FRAMEWDRK    AND 
GENERALIZED    BUDGET   EQOATIDNS 

The  purpose  of  this  appendix  is  to  present  the  formula- 
tion of  the  storm  budget  volume  coordinate  system  and  the 
generalized  budget  equations  which  can  be  applied  to  any 
desired      atmospheric   property.  The   information     contained 

herein  is  extracted  from  Wash  (1978)  as  adopted  from  Johnson 
and    Downey    (1975   a). 

The  storm  budget  volume  in  spherical  coordinates  is 
depicted  in  Fig.  4  1.  The  follDwing  notation  is  relevant  to 
this    figure: 


a 


Azimuthal   coordinate 
Angular  radial  coordinate 


Position   vector  from    earth's 
center 


Position  vector   from    storm   center 
to   any    point    in   budget    volume 


a,l,k  Orthognal   unit   vectors    along 

~  «-  "-  meridional,    azimuthal   and 

vertical  directions   of    spherical 
coordinates    for   storm 


Mind   velocity   relative    to   earth 


1  13 


Horizontal  velocity  of    reference 
axis 


Horizontal  velocity  of   any  point 
on  the   volume  boundary. 


Table  71  contains  a  listing  of  the  generalized  budget 
equations.  The  following  is  an  explanation  of  the  notation 
used   in   these  equations: 

a  The  vertical    coordinate   system 

V  Volume 

T  Top  boundary 

3  Bottom  boundary 


(U-W)        Horizontal  wind  component  relative 
0  to  budget  volume  normal  to  lateral 

boundary 


(drr/dt)      Vertical  velocity  divergence  into 
a  layer  within  the  volume 


rsin/2        Radius  of  budget  voluie 


p  Surface  pressure 


Specific  budget  property  (per  unit  mass) 
Vertical  velocity  (d<7/dt) 
Acceleration  due  to  gravity 
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Earth's 
surface 


Figure  4 1. 


Quasi- Langrangian  storm  budget  volume  coordinate 
system. 
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TABLE    V 

GENERALIZED    QDASI-LAGRANGI AN    BUDGET    EQUATIONS    IN    SIGMA 

COORDINATES 
(after  Wash,    1978) 


The   budget   volume  integral   is   defined  as 
a       B„     2tt 
F  «  fr     sinBdadB(-da). 

aB     0       0 
The   most    general    form,  of    the  budget    equation  is 

dF 

^-     »     LT   (F)        +       VT(F)        +       S(F) 

at 

where  the   lateral  transport  term  is 
LT(F)     -J  — -  f(U-W)gr  sin  Bda(-da)    |BB  , 

the  vertical  transport   term  is 

a-    h    2* 

p+  • 

r-  (fa)  sinSdadB  (-da) , 


VT(F)      -J  /  /  J^*     Z_ 

aB        0       J0  S        3a 

and  the   source    (sink)    term   is 

a      ^B      2tt 

S(F)     "J     J    J  T  |fr2sInMod6(-da)- 

°B       o         0 

The    vertical   transport  tern   is    zero 
when   the   limits   of " i ntegration    are 
<J  -    1    and     (7=0,  since   there   can   be 
no   transport    through   the   top   and 
bottom   boundaries  of   the  budget   volume. 
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